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A systematic study of the properties of particle and charge correlation functions in 
the two-dimensional Coulomb gas confined to a one-dimensional domain is under- 
taken. Two versions of this system are considered: one in which the positive and 
negative charges are constrained to alternate in sign along the line, and the other 
where there is no charge ordering constraint. Both systems undergo a zero-density 
Kosterlitz-Thouless type transition as the dimensionless coupling T := q 2 /kT is var- 
ied through r = 2. In the charge ordered system we use a perturbation technique 
to establish an 0(1/ r ) decay of the two-body correlations in the high temperature 
limit. For r — > 2 + , the low-fugacity expansion of the asymptotic charge-charge 
correlation can be resummed to all orders in the fugacity. The resummation leads 
to the Kosterlitz renormalization equations. In the system without charge ordering 
the two-body correlations exhibit an 0(l/r 2 ) decay in the high temperature limit, 
with a universal amplitude for the charge-charge correlation which is associated with 
the state being conductive. Low fugacity expansions establish an 0(l/r r ) decay of 
the two body correlations for 2 < L < 4, and an 0(l/r 4 ) decay for L > 4. For 
both systems we derive sum rules which relate the long wavelength behaviour of the 
Fourier transform of the charge correlations to the dipole carried by the screening 
cloud surrounding two opposite internal charges. These sum rules are checked for 
specific solvable models. Our predictions for the Kosterlitz-Thouless transition and 
the large-distance behaviours of the correlations should be valid at low densities. 
At higher densities, both systems might undergo a first-order liquid-gas transition 
analogous to the two-dimensional case. 
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l.INTRODUCTION 



The natural domain for a system of charges interacting via the two-dimensional (loga- 
rithmic) Coulomb potential is a plane. Notwithstanding this fact, there is still interest in 
studying the statistical mechanics of log-potential Coulomb systems confined to a line (we 
will refer to such systems as log-gases). Two-component log-gases have attracted interest 
because of some equivalences with models in solid state physics (Kondo problem^ 1 ) and 
quantum Brownian motion problem^ 2 ) ), and the fact that for low density a Kosterlitz- 
Thouless-type pairing transition takes place as the temperature is lowered below the 
critical coupling T = 2^' 3 \ 

In this paper we will study properties of correlations in the two-component log-gas 
with oppositely signed charges of strength q, and a variant of this system in which the 
positive and negative charges are confined to alternate in position along the line. For 
couplings T := q 2 /kT > 1 we also assume that the particles are at the centre of hard 
rods of length a , which prevents short distance collapse. So as to put our work in 
context, let us briefly review known properties of the critical behaviour and correlations 
in two-component log-gas systems, and contrast this with what is known about the two- 
component two-dimensional Coulomb gas (2dCG) where appropriate. 

1.1 The charge ordered system 

The system with charge ordering was first studied in the context of its application to the 
Kondo problem^ 1 ). In this seminal work, Anderson et al. transformed the grand partition 
function of the system with a hard rod length a + da, fugacity ( and coupling T into 
the grand partition function of the system with hard rod length a and modified ( and T, 
thereby deriving a pair of coupled renormalization equations. The transformation, which 
is approximate, requires ( and T— 2 to be small and is thus applicable in the neighbourhood 
of the Kosterlitz-Thouless transition. Remarkably, applying this procedure to study the 
2dCG , Kosterlitz^ 4 ) found essentially the same equations. Furthermore, for the 2dCG it 
is well known (see e.g. ref.[5 ]) that a renormalization procedure can be applied to study 
the charge-charge correlation and a certain length dependent dielectric constant. Again, 
by an appropriate choice of variables, the resulting equations are precisely those found by 
Anderson et al.. 

One immediate prediction from the renormalization equations is the dependence of 
the critical coupling (r c ) on fugacity: thus for the charge ordered system 

r c -2 = 2 5/2 c (l.i) 

valid to leading order in T c — 2 and (. This phase boundary separates a low-temperature 
dipole phase in which the positive and negative charges are paired, from a high tempera- 
ture phase in which the positive and negative charges are dissociated. More quantitatively, 
from the mapping to the Kondo problem, Schotte and Schotte^ have argued that for the 
finite charge ordered system of length L with periodic boundary conditions 

//V^ r /r>\2\ f C (F)L for T < 2 

((E^-W^U^ forr>2 ^ 

This behaviour, in which the boundary conditions play an essential role, was confirmed 
by Monte Carlo simulation^ 6 ) . 
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In contrast to the phase indicator (1.2), the phases of the 2dCG can be distinguished 
by the dielectric constant e, which is finite in the dipole phase and infinite in the high 
temperature (conductive) phase. We recall that if the system has dielectric constant 
e, then a fraction 1 — 1/e of an infinitesimal external charge density will be screened. 
Furthermore, e can be expressed in terms of the particle correlations by the formulas*- 7 -* 

I = 1 + ^ f drr 2 C(r) (1.3a) 



and 



(3q 
2p 



J drF(r)p- + (r)p- + (r) (1.36) 



— * 

where C{r) denotes the charge-charge correlation, F(r) denotes the force corresponding to 
the two body potential between two positive charges of unit strength (which is assumed to 
be smoothly regularized at the origin), p~+(r) denotes the distribution function between 
two opposite charges and p_+(r) denotes the dipole moment of the charge distribution 
induced by two fixed opposite charges separated a distance r. 

In log-gas systems it can be argued that a dielectric phase has e = 1 ( see e.g. ref. [8]). 
A linear response argument then gives 

C(r) ~ o (-^j as r ^ oo (1.4a) 

On the other hand, in a conductive phase the charge-charge correlation must exhibit the 
universal asymptotic behaviour*^ 

C(r) ~ -pgp d.46) 

In the sense of screening an infinitesimal external charge, the charge ordered system is 
expected to always exhibit a dielectric phase (except possibly at the phase boundary) ®. 
Hence, from (1.4a), the decay of C(r) should be faster than 1/r 2 in both the low tempera- 
ture dipole and high temperature phase. When the positions of the particles are restricted 
to the sites of a lattice (this mimics a hard-core regularization, the isotherm T = 1 is ex- 
actly solvable^ 9 ) . Both the truncated two-particle distribution functions between like and 
opposite charges, and thus C(r), exhibit an asymptotic 0(l/r 4 ) decay in accordance with 
this prediction. 

1.2 No charge ordering 

Application of the rescaling method of Anderson et. al. to the two-component log-gas 
without charge ordering leads to the conclusion® that the critical coupling is T = 2 
independent of the fugacity, in contrast to the behaviour (1.1). No other information is 
obtained. 

The charge-charge correlation function in the low temperature dipole phase must again 
exhibit the behaviour (1.4a). This has been explicitly verified for the solvable isotherm 
T = 4, where 0(l/r 4 ) decay was found^ 10 ). The high temperature phase is expected 
to be conductive and thus C(r) should obey the sum rule (1.4b). The phase boundary 
T = 2, when the positions of the particles are restricted to the sites of a lattice, is a 
solvable isotherm* 10 ' 11 ) and it is found that C(r) has a 0(l/r 2 ) decay, but with a density 
dependent amplitude. 
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1.3 Paper outline 

From the above brief review, it is clear that there are many gaps in our knowledge of the 
critical properties and the behaviour of correlations in two-component log-gas systems. 
To improve on this situation, we will make a fairly systematic study of the correlation 
functions in the high and low temperature phases as well as in the vicinity of the zero- 
density critical point. The charge ordered system is considered in section 2, while the 
system without charge ordering is considered in section 3. Concluding remarks are made 
in section 4. 

In subsection 2.1 the two-particle correlations in the high temperature phase of the 
charge ordered system are analysed using a perturbative approach. The correlations in 
the scaling region are analysed in subsection 2.2 by applying the low density resummation 
method of Alastuey and Cornu^ and in subsection 2.3 the correlations in the low tem- 
perature phase away from criticality are considered. In subsection 2.4 the second BGY 
equation is transformed into Fourier space and a sum rule analogous to (1.3b) is obtained. 

In subsection 3.1, the two-particle correlations in the conductive phase of the system 
without charge ordering are studied using a linear response argument and macroscopic 
electrostatics, and a high temperature resummation technique. The asymptotic behaviour 
of the correlations in the dielectric phase are determined in subsection 3.2 by studying 
the low fugacity expansions at 0(( 4 ), and an analysis of the second BGY equation similar 
to that given in subsection 2.3 for the charge ordered model is made in subsection 3.3. 
From the latter analysis the behaviour of the three and four particle correlations in the 
dielectric phase is deduced. 

In Section 4 we summarize our results with emphasis placed on the mechanisms behind 
the contrasting behaviours of the correlations in the log-gas with and without charge 
ordering. 

2. THE CHARGE ORDERED LOG-GAS 

2.1 Decay of correlations in the high temperature phase 

Along the high temperature solvable isotherm T = 1 it has been shown by explicit 
calculation^ 9 ) that the charge-charge correlation decays as 0(l/r 4 ), and thus by (1.4a) the 
system does not screen an infinitesimal external charge. Indeed, because of the charge 
ordering constraint, it was argued^ 9 ** that this latter property is a general feature of the 
system for all temperatures. In particular, the high temperature phase is not conductive, 
which suggests that the conventional methods of analysing the high temperature phase 
in Coulomb systems (Debye-Huckel type theories etc.) are not applicable. In fact the 
familiar Abe-Meeron* 12 ) diagrammatics, with Debye-like mean-field potentials resulting 
from chain resummations, cannot be applied here because the constraint of charge order- 
ing is equivalent to introducing a many-body potential between the charges, whereas the 
diagrammatics hold for systems with two-body forces only. We use instead a perturba- 
tive method which relies on the special screening properties of the T = 0, constrained 
two-species perfect gas reference system. 

2.1.1 Correlations at T = 

For T < 1 in general, and T = in particular, there is no need to regularize the short 
range singularity of the logarithmic potential as the corresponding Boltzmann factors are 
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integrable. Thus at T = the log-gas system can be considered as a two-species perfect 
gas of point particles, constrained so that the two-species alternate in position along the 
line. Let us suppose the first particle belongs to species + while the second particle 
belongs to species — (these boundary conditions will result in a symmetry breaking: the 
p + _ and p_ + correlation functions will not in general be equal), and there are N particles 
of each species with coordinates x±, . . . , x^ (+ species) and yi, ■ ■ ■ ,Un (— species) where 
Xj, %]j E [-L/2, L/2](j = 1, . . . , N). The canonical partition function Z 2N for this system 
is defined as 

Z2N-=(^7-.) / dx l ...dy N (2.1a) 



where * denotes the interlacing constraint on the interval [-L/2, L/2]. It is easily evalu- 
ated to give 

z - = my 1 *" (m 

The two particle distribution p ++ (x 1 ,x 2 ) is given by 

. , N(N - 1) /■*(si>*2) 

p++{x 1 ,x 2 ) := — jz / dx 3 ...dy N (2.2a) 

where 

Q 2N = (N\) 2 Z 2N (2.26) 

and *(xi,x 2 ) denotes the interlacing constraint, given that there are particles of species 
+ at x\ and x 2 . The distribution p_ + (xi,yi) is defined similarly. 

To calculate p ++ (xi,x 2 ), we note that in the interval [—L/2,xi] there must be equal 
numbers of particles of each species, while in the intervals [x\, x 2 ] and [x 2 , L/2] there must 
be one more particle of species — than particle of species +. Let the number of particles of 
species + in each interval be M a , Mb and M c respectively (note that M a +Mb+M c = N—2). 
Then we have 

, _ (AH) 2 ^ 2 (xi + L/2) 2 ^(a; 2 - Xl f M ^\L/2 - x 2 ) 2M ° +1 ( . 
P++[Xl ' X2) ~ Z 2N Ma ^ (2M a )!(2M fc + l)!(2M c + l)! ^ 

M a +M b +M c =N-2 

Substituting (2.1b), this can be rewritten as 

p++(x 1 , x 2 ) = 2N ^ 2N ^ >S , 2jv-2(^i + L/2, x 2 - Xi, L/2 - x 2 ) (2.4a) 

where 

p/2-i I 

O I U \ ._ V 2l 2M af 2M b +l 2M C +1 icy Ah \ 

S * aM - M ^ Mc (2M a )!(2M b+ l)!(2M c + l)! a b ° {2M) 

M a +M b +M c =p/2-l 

Using the generalized binomial expansion 

p v \ 

(a + b + cf= Y ,\ a pi F 2 c Pi (2.5) 

P1+P2+P3=P 
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it is straightforward to derive the summation formula 

S p (a, b, c) = i [(a + b + cf + (-a + b + c) p -(a-b + cf -(-a-b + cf] (2.6) 

valid for p even. With this result the thermodynamic limit in (2.4a) can be taken imme- 
diately to give 

P ++ (0,a;)=p 2 [l-e- 4 ^] (2.7) 

A similar calculation shows 

p + _(0,x) = p 2 [l +e- 4 ^l] (2.8) 
or alternatively this result could be deduced from the requirement 

^( P++ (0,x)+ P+ 40,x))=p 2 (2.9) 

which follows since the combination of two particle distributions on the l.h.s. gives the 
two-particle distribution of the (unconstrained) perfect gas. 

The crucial feature of these distributions is that they exhibit perfect screening of an 
internal 'charge'. Thus 

/oo 
dx (p ++ (0, x) - p + _(0, x)) = -p (2.10) 
-oo 

This property allows the correlations in the high temperature limit to be studied by a 
perturbation expansion about the T = constrained perfect gas reference system. 

2.1.2 Perturbation about T = 

In any two component log-gas system, the Coulomb interaction energy U can be written 



U = \j dx j dx'[Q(x)Q(x')} nc v c (\x-x'\) (2.11a) 



where 

v c (\x-x'\) := -\og\x-x'\ (2.116) 

N 

Q(x) := g£ ~ Xj) - S(x ~ Vi)) (2-Hc) 

3=1 

with the positive (negative) species at Xj(yj) and 

[Q(x)Q(x')] nc 

denotes that the product Q(x)Q(x f ) is formed with products over coincident points ex- 
cluded. In the charge ordered two-component system the truncated two particle distribu- 
tion between like charges is given by 

t rn , _ (iN + (0)N + (x a )L c e~n (N + (0)e-^)o(N + (x a )e~no (9 

where the subscript indicates that the averages are taken with respect to the constrained 
perfect gas reference system and 

N 

N + (x a ) = J2$( x ~ x j) 

3=1 
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Next we expand the exponentials in (2.12) to leading order in (3. We obtain 

R roo roo 

p T ++ (0,x a ) ~ -£ / / dxdx'(N + (0)N + (x a )[Q(x)Q(x')l c )W^-^\) (2-13) 

Z J —oo J — oo 

where the truncation is defined with respect to the three quantities 0,x a and (x, x')S 13 ^ 
When expressed in terms of the fully truncated Ursell functions there are two contribu- 
tions. One involves only two body Ursell functions while the other involves the four body 
Ursell function. Since the Ursell functions in the reference system decay exponentially, 
the leading contribution comes from the term involving the two body functions. Hence 
(2.13) can be rewritten as 

/oo roo 
/ dxdx' (N + (0)Q(x)} (N + (x a )Q(x')} v c (\x - x'\) (2.14) 
-oo J —oo 

The two body averages are given by 

(N + (0)Q(x)) = q (pf + (x) - p {0 l(x) + p S(x)) = -2qp 2 e- A ^ + qp 5(x) (2.15a) 

and 

(Q(x')N + (x a )) = (Q(x' - x a )N + (0)) (2.156) 
To calculate the leading large-x a behaviour of (2.14) we make the expansion 

oo / / _ _ \n On 

v c {\x -x'\)= v c {\x a \) + E [ ~ ^-^c(k|)|_ o (2.16) 

Due to the perfect screening property (2.10) of the reference system, and the fact the 
distribution functions in the reference system are even in their argument, we see that the 
first non-zero term which results from substituting (2.16) in (2.14) is n = 4. Noting 

d 4 „ , x <9 4 , . 3! 



i u c(N) = -^-7log|a; 



x 4 



dx 4 dx 4 
we therefore have that to first order in T and for large-x a 

~tS& <2 - 17) 

The leading asymptotics of p+_(0, x a ) can be computed in a similar way. In fact, since 
in the reference system 

pl_(0,x) =-pl + (0,x) 
the only difference in the calculation is a minus sign so we obtain 

pU 0, x) ~- pM ~^ (2,8) 

Although (2.17) and (2.18) are valid to first order in T only, we expect the leading 
0(l/a; 4 ) decay to persist in the high temperature phase for T < 1 at least, as the exact 
result at T = l^ 10 ) exhibits this behaviour. Indeed we expect all terms in the T-expansion 
of the truncated two body distributions to decay as 1/x 4 ; this would follow from a di- 
agrammatical analysis similar to the one used by Alastuey and Martin^ 13 ) . The exact 
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result also exhibits the property (2.18) relating the leading asymptotics of p+_(0,:r) and 
P++(0, x), which suggests this may also be a general property of the high temperature 
phase for T < 1 at least. 

2.2 Correlations in the scaling region of the low temperature dipole phase 

For the 2dCG it has been proved^ 14 ** that all the coefficients in the low fugacity expansion 
of the pressure and the correlation functions are convergent in the dipole phase (r > 
4). We expect this property to remain true in the low temperature dipole phase of 
the log-gas systems (r > 2). With this assumption, we can consider the asymptotic 
large-r behaviour of each such coefficient in the charge-charge correlation C(r) for T > 
2. In particular we can consider the behaviour for r — > 2 + , as the Kosterlitz-Thouless 
transition is approached from the dipole phase. Motivated by the analogy between critical 
features of the alternating model and the 2dCG (in particular the occurrence of the same 
renormalization equations), we do this by closely following the strategy of Alastuey and 
Cornu^ in their analysis of C(r) in the two-dimensional system. 

For this purpose we introduce the asymptotic charge density C&(r). This is defined 
as the terms in the asymptotic expansion of C(r), which when replacing C(r') in 

A:=l + — / dr'r'C(r') (2.19) 

give the correct leading order singular behaviour of A for T — > 2 + at each order in (. Note 
from (1.3a) that A is the analogue of 1/e in the 2dCG. Furthermore, in the low fugacity 
limit A — 1 is proportional to the mean distance of separation between neighbouring 
charges* 15 ^. Our objective is to calculate Ca(^) and A in the scaling region of the low 
temperature dipole phase. 

2.2.1 Low fugacity expansions 

The low fugacity expansions of the correlations can be obtained from the low fugacity 
expansion of the logarithm of the grand partition function, with each positive (negative) 
charge given a position dependent fugacity (a(xj)((b(yj)). The grand partition function 
is given by 

oo 

E[a,b]= Y.C N Z 2N [a,b] (2.20) 
n=o 

where Z 2 ^[a, b] denotes the partition function for equal number N of positive and negative 
charges. The truncated distributions can be calculated from (2.20) according to the 
formula 

...,r ni ;si,...,s n J = — — — — — - log S[a,o] 

(2.21) 

Now from (2.20) 

log~M = C 2 Z 2 [a,b] +C 4 (z A [a,b] - l -{Z 2 [a,b]f) + 0(C 6 ) (2.22) 

Using (2.22) in (2.21), assuming that the left most particle always has a positive charge, 
and writing the partition functions explicitly we obtain for the two-particle correlations 



pl + {r)=C 



rr—a roo 

/ dy 1 / dy : 

Jit Jr+a 



r(V2 - Vi) 



r+v \(y2 - r)y 2 (r - yi)yi 




+ 0(C 6 ) (2.23a) 
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dx 



+2 



2a 



dy 



(r — y)x 



r(x — y)(r — x)y 



°° dy 



roc fy~°~ 

I dy / 

r+2a J a+r 



dx 



x(y — r) 



(y — x)y(x — r)r 

| / rr+2a ry — (7 roc rr 

f ( / dy dx + dy dx 

T \J —o J —oo Jr+a J —oo 



(y - x) r 



(y - x)r) 

+ 0(c e 



r) = C 4 



roo r — cr 

/ dy 

Jr+cr J —oo 



dx 



(r - x)y 



r(y — x)(y — r)x 



00 dy 



(2.236) 



+ 0(C 6 ) (2.23c) 



where it is assumed r > 2a in (2.23a) and r > a in (2.23b) and (2.23c) (for r less than these 
values the respective full distribution functions vanish). Also, the first double integral in 
(2.23b) is to be omitted if r < 2a. To the same order, the low fugacity expansion of the 
charge-charge correlation follows immediately from the formula 



C(r) = 2q 2 p6(r) + g 2 [2p^(r) - p^(r) - p__ + (r)\ 



(2.24) 



2.2.2 Evaluation of C 



(2), 



r) and C^(r) 



Let us denote the term proportional to C 2 - 7 in the low fugacity expansions of C(r), C&{r) 
and A by C^\r),C^ ]) {r) and A (2i) respectively. From (2.24) and (2.23) we have, for 
r > a, 

C^( r ) = q 2 ^ (2.25) 

Substituting in (2.19) gives 

A m = i + 4r f V ~ F (2.26) 



r -2 

We note that A^ 2 ^ is singular as T — > 2 + and furthermore to leading order is independent 
of a. Both features are true of A^ 2 ™) in general. The latter feature implies that only the 
large-r portion of C^ 2n \r) contributes to the leading order singular behaviour of A^ 2n \ 
and thus C^ n \r) consists of terms in the asymptotic expansion of C^ 2n \r). With n = 1 
there is only one term in the asymptotic expansion, which is C^(r) itself, so trivially 
C { l\r) = CW{r). 

The analysis of the large-r behaviour of the integrals in (2.23), which from (2.24) form 
C^(r), is done in appendix A. There it is deduced 



;r - 2)- 



r / 



2-r 



1 



r -2 



2-r 



logr 



(2.27) 



and the integral representation 

. . 1 r2r r2x—r 

C ( l\r) = -q 2 C- dx dyS {0 . 

T Jr+a Jx+a 

AY 

S(o,r)(x,y) ■= 



r) 



[x,y) 



with 



,x 



r)(y-x) 



r-i 



(2.28a) 



(2.286) 
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is also derived. The r.h.s. of (2.27) has the remarkable property of having an identical 
structure to the corresponding expression in the low fugacity expansion of the charge- 
charge correlation in the 2dCG [ref. 7,eq.(4.13)]. Furthermore, the integral expression 
(2.28) can be interpreted as the contribution of a mobile positive-negative dipole (positive 
charge at x, negative charge at y) pair partially screening (note that in all phases the 
external charges are not screened at all since e = 1) the root dipole of separation r 
(positive charge at 0, negative charge at r). The distance of separation \y — x\ between 
charges within the mobile dipole is constrained to be less than the closest distance between 
this dipole and the root charges at or r, and this closest distance is to be no greater 
than r. The factor of 4 comes from the 4 equivalent ways of arranging the mobile dipole 
about the root charges (the screening dipole may lie close to or r, and inside or outside 
the root dipole). This "nested" dipole interpretation of C^\r) is analogous to that found 
in ref. [7] for the quantity C e (4) (r) in the 2dCG. 

2.2.3 Nested dipole chain hypothesis 

Although we have calculated C^\r), it may seem a formidable task to do likewise for 
Cf J '(r), j > 3. Indeed, an analysis similar to appendix A does not appear to be feasible. 
Instead, having identified an analogy between C^\r) and C^ 4 \r) of the 2dCG, we adapt 
the method given in ref. [7] to calculate Cf j \r). 

The basis of the method of ref. [7] is the hypothesis, which has its origin in the work of 
Kosterlitz and Thouless^ 16 \ that the configurations contributing to C^\r) are all nested 
chains of dipoles, with the fixed dipole (positive charge at and negative charge at r) 
the largest. The screening operator (2.28b) acts between dipoles connected in a chain. 
Furthermore, by including the factor of 4 in (2.28b) these chains can all be ordered to 
the right of the fixed negative charge at r. For example, at 0(( 6 ), there are two distinct 
chains as given in Fig. 1. 

The contributions to C^^r) from these chains are 



r2r r2x 2 -r r"i'ij2-^2 f%x\—y2 

/ dx 2 / dy 2 S^ r) (x 2 ,y 2 ) / dx 1 / dy^^y^x^yi) 

Jr+a JX2+(T Jy2+a Jxj+a 



and 



r2r r2x2—r 1 r r'Zr r'Axi~r 

/ dx 2 / dy 2 S^ r) (x 2 ,y 2 ) / dx x \ dy x S^ r) {x x ,y x ) 



r2r f2x\—r 



Furthermore, the first of these contributions needs to be weighted by a factor of 2 to 
account for the relabelling degeneracy, and this linear combination needs to be multiplied 
by a factor of — g 2 ( 6 /(4!r r ). 

At general order, the nested chain hypothesis gives [ref.7 ,eqs.(4.26)-(4.28)] 

2f2 A(2n-2) 

CfV) = -7r-(^32)!^ 2n - 2) (r) (2.29a) 

where 

o(2n-2),x ^ (W-I)! ^ (71- 1 -p)\ 



qi+...+q p =n-l-p 



with 



P 2 T P^X — 7" 

I 2q (r) = dx dyS {0>r) (x,y)S < £ q \x - y) (2.29c) 

Jr+a J x+a 
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From the structure (2.29) it is shown in ref. [7] that the series 



Ca(t) = £ c A 2 



(2n), 



(2.30) 



n=l 



can be summed whatever the form of S(p tr )(x, y), with the result 



Ca(t) 



q\ 2 

— — exp 



2r 



dx 



Jr+u 



2x—r 



x+a 



dyS(o, r ) (x, y) (y -x) T C A (y- x) 



(2.31) 



Since S(o, r ){x,y) is explicitly given by (2.27b), we see that it is possible to simplify the 
double integral in (2.31) by an integration by parts in the y-variable. This gives 



q 2 C 2 
— — exp 



Ca(t) 



rr rr 

-4(3 log r / dxxC A (x) + 4/3 / dx \ogxxC A {x) 

J cr J a 



(2.32) 



This expression has the same structure with that of C e {r) for the 2dCG [ref. 7,eq. (4.32)]. 
Introducing the length dependent version of (2.19) by 



AT r 

A( r ) := 1+ / dr'r'C A (r') 
q 2 Jo 



(2.33) 



which gives the contribution to A from all particles of separation at most r and is anal- 
ogous to the length dependent dielectric constant defined by Kosterlitz and Thouless^ 16 - 1 , 
we can obtain from (2.32) and (2.33) the coupled differential equations 



ctlog(r) q 2 



and 



d 



d\og(r 

Making the change of variables 



-C A (r) = -rC A (r)A(r) 



(2.34a) 



(2.346) 



y(r) = - (-rC A (r)) 1/2 r, t(r) = 1 - ^A(r) 
q Z 



I = logr 



the differential equations (2.34) read 



d_ 
Jl 

d 
Jl 



y(i) = y(i)t(i) 



t(i) = {y{i)Y 



(2.35a) 



(2.356) 



These equations, with different meanings to y{l) and £(/), are precisely those obtained by 
Anderson et al.^. 

It is simple to solve the system (2.35) for £(/) in terms of y(l), where the constant of 
integration can be determined by putting r = a in the integral equation. We find, using 
r instead of /, 



-^C A (r) = (l- l -A(r)y + 4IXV 



2 2-r 



Y \ 2 
1_ 2, 



(2.36) 
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Now, assuming T > 2, from the low fugacity expansion Ca(V) is o(l/r 2 ) as r — > oo. 
Recalling the definition (2.19) and taking the limit r — > oo in (2.36) we thus have 



T-2 
1 _ 



|-(-K)V- n1 
(r-2) 2 



- 1 (2.37) 



Notice that the radius of convergence of the resummed function in (2.37) gives the phase 
boundary (1.1). This expression for A has an identical structure to that of 1/e in the 
scaling region of the dipole phase of the 2dCG ( - 7 ^. 

The resummation of C^(r) can be performed by analysing the differential equations 
(2.34). This has been done in ref. 7 using these equations as they apply to the correlations 
of the 2dCG. From the workings of ref. 7 [eqs.(4.41)-(4.50)] we find 



Ca(0 = -q 2 C 

where 



/<7\ rA * /a\ rA+w(rA-2)\ 

^o(-) +E^( 7 ) ) (™°) 



A„ = exp|-r£ °*[A(()-A]J (2.386) 
and the A^r can be expressed in terms of A and the function 

2(4C) 2 /(r-2) 2 



[l-2(4C) 2 /(r-2) 2 ]V2 



(2.39) 



Thus the coupling T is renormalized as the density dependent quantity TA, and TA is 
the power involved in the expansion (2.38) of C&{r). 

2.3 Correlations in the low temperature dipole phase away from criticality 

The analysis of Appendix A identifies the leading large-r expansion of p+ + (r),p^_(r) and 
p T + (r) at 0(( 4 ) throughout the low temperature dipole phase T > 2. Thus from the 
expansions (A6) we have 

P lf(r) ~ p?_ 4) (r) ~ P T _f{r) ~ (2.40) 

In the case of p++^(r) this behaviour results from the contribution of the configuration A 
of figure 2 in Appendix A. Similarly, configurations in which each root particle is paired 
with a mobile particle to form a neutral cluster give the leading contribution to p+^(r) 
and pT^(r). These configurations have the same statistical weight at large distances, 
which implies the same behaviour (2.40) for p*+- and p 7 ^ . 

To calculate the large-r behaviour of p++ n ^(r) etc. for n > 2 we need to hypothesize 
the generalization of configuration A which is dominant in the respective integral formulas 
for r — > oo. This generalization is to have all particles belonging to a neutral cluster about 
and including one or other of the root particles, with the interparticle spacing within a 
cluster small in comparison to r. The potential V 2 n between clusters can be written as 

V 2N = W + W r + U 0r (2.41) 
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where Wq and W r denote the electrostatic energies of the individual clusters about and 
including the root particles at and r respectively, and Uq t denotes the mutual interaction. 
For large-r, U 0r is to leading order a dipole-dipole potential: 



d d 



dx dr 



PoPr 



Uor ~ POTT-Pr^ log |x - A _ = (2.42) 



where p and p r are the dipoles of the clusters about and including the root particles. 
The situation for n > 2 is thus completely analogous to the case n = 2 (recall the text 
below (A2)). Insertion in the corresponding integral formulas, and use of the expansion 

e-PUor „ i_p U()r + _ m (2.43) 

thus leads to the results 



P++ ( r ) ~ P+- ( r ) ~ P-+ ( r ) ~ 2 ( 2 - 44 ) 

where the coefficient a 2n (r) diverges as T — > 2. This 0(l/r 2 ) behaviour can be under- 
stood as resulting from the dipole-dipole interaction between the neutral clusters at large 
separation. 

The configurations giving the leading large-r behaviour of higher order correlations 
will again be neutral clusters and, analogous to the case n — 2, the potential between far 
away neutral clusters will again be to leading order due to the dipole-dipole interaction 
and thus 0(l/r 2 ). Hence we expect 

^...-<*.....,*«;».,...,o ~ n^n^ni^-wi' (2 ' 45) 

as Xjk := Xk — Xj, y^w '■= Uk' — Uj', \%j — Vk'\ where we have omitted the amplitude. 

This asymptotic behaviour requires at least one mobile particle about each root particle, 
and thus does not necessarily apply for n < n\ + n 2 . 



T(2n)/ x T(2n)/ x ^ C W «2n(r) 



2.4 A sum rule from the BGY equation for C(x) 

It has been shown in ref. 7 that the BGY equation for the charge-charge correlation C{r) 
in the 2dCG can be used to derive a sum rule which expresses the dielectric constant 
in terms of the dipole moment p+~(r) of the screening cloud surrounding two internal 
charges of opposite sign. In this subsection we will apply an analogous analysis to C(x) 
(we use x to denote any position coordinate, positive or negative and r to denote a non- 
negative quantity) for the charge ordered system which leads to a sum rule involving the 
distribution functions between charges of opposite sign and the dipole moments p + -(x) 
and 

The BGY equation for C(r) in the 2dCG has been derived from the BGY equations 
for the two-particle distributions. The BGY equation for C(x) in the two-component 
log-gas without charge ordering can immediately be read off from the 2dCG result (ref. 7, 
eq.(5.10)). However, the resulting equation is not applicable to the charge ordered system 
for two reasons: (i) the charge ordered constraint introduces extra terms in the BGY 
equation for the two-particle distributions and thus for C(x), and (ii) the symmetry 
p + _(x) = p_ + (x) assumed in the derivation is not valid in the low temperature phase of 
the charge ordered system. 
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To derive the BGY equations for C(x) in the charge ordered system let us then re- 
consider the BGY equations for the two-particle distribution functions. We suppose the 
logarithmic potential is regularized by hard cores of length a symmetrically placed about 
each particle. We know that without charge ordering, and with smoothly regularized 
potentials 

= s 1 s 2 v{x) 

the BGY equation for the distribution pJ 1S2 (x) can be written in the form (ref. 7, eqs. 
(5.8a), (5.8b)) 

q^P^is2 faia) = SiS2TF(x 12 )p^ lS2 (x 12 ) + s 2 f3qp 2 J dx 3 F(x 32 )Q Sl (%i \x 3 ) 

/oo 
dx 3 F(x 32 )[pJ lS2+ (x 1 , x 2 , x 3 ) - p T SiS2 _(x u x 2 ,x 3 )} (2.46a) 
-oo 

where 

d 

F(x) := -—v(x), x ab := x b - x a (2.466) 

and Q sl ... Sn ( 

) denotes the total charge density induced at x given that there 
are charges si, . . . , s n fixed at x±, . . . , x n \ 

f-\ i I \ [Psi...s n +{%li ■ ■ ■ i x ni Psi...s n — ■ ■ ■ i x ni %)\ 



Psi...s„(%li ■ ■ ■ i %n) 
n 

+qJ2siS(x-Xi) (2.46c) 
i=i 

By considering the definition of p SlS2 (x) in the canonical ensemble, it is easy to see that 
the only modifications needed to this equation due to the charge ordering constraint and 
the hard cores are that 

p slS2 (- S2 )(x 1 ,x 2 ,x 2 - a) - p SlS2 (- S2 )(x 1 ,x 2 ,x 2 + a) (2.47) 

must be added to the r.h.s., and the smoothed potential v(x) replaced by v c (x) (recall 
(2.11b)). The definition of Q si (xi\x 3 ) remains that given by (2.46c) and the integration 
over x 3 remains extended over the whole real line. This is a consistent prescription as 
the r.h.s. of (2.46a) can then be expressed entirely in terms of full distribution functions 
which vanish for configurations with — Xj\ < a. However, the equation (2.46a) modified 
by the hard core terms (2.47) is now valid for \x± 2 \ > a (notice that the l.h.s. of (2.46a) 
diverges at \x± 2 \ = a) 

Adding (2.47) to the r.h.s. of (2.46a), multiplying both sides by q 2 S\S 2) summing over 
Si, s 2 — ± and use of the definitions (2.24) and (2.46c) gives that the BGY equation for 
C (x) in the charge ordered system with hard cores is 



d 



(C(x 12 ) - 2q 2 p5(x 12 )) = 2Tp f dx 3 F c (x 32 )C(x 13 ) + Tq f dx 3 F c (x 32 )R(x 1 ,x 2 ,x 3 ) 



dx 2 

+qA( Xl ,x 2 ) (2.48a) 

where 

R(x 1 ,x 2 ,x 3 )= J2 Q s i s 3pJ 1 s 2 s, i ( Xl r x 2,x 3 ) + qJ2p^ 1 s 2 ( x u x 2)S(x 3 -x 1 ) (2.486) 

S1,S2,S3=± S1S2 
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= s 3 p S2S3 (x 2 ,x 3 )Q S2S3 (x 2 ,x 3 \x 1 ) - — C(x 31 ) (2.486') 

S 2 ,S 3 =± 1 

(the symbol = denotes that (2.48b') gives the same value as (2.48b) when substituted in 
(2.48a); see the sentence below (2.48d)) and 

A(xi, x 2 ) = B(xi ) x 2) x 2 — a) — B(x\ ) x 2) x 2 + a) (2.48c) 

with 

B(x 1 ,x 2 ,x) = p + _(rr 2 , x)Q + ^(x 2 , x\x±) - p_ + (rr 2 , x)Q -+(x 2 , x\xi) 

-q(5(xi -x 2 ) -8(xi - xf)(p+-(x 2 ,x) + p^+(x 2 ,x)^ (2AM) 
In deriving (2.48) we have used the equation 



/oo 
dxF c (x) S 2 S 3 p 8283 (0,x) = 

-oo „_ L 



S2,S 3 =± 

which is a consequence of the first factor in the integrand being odd while the second 
factor is even. 

In preparation for analysing (2.48a) we note that on the l.h.s we can make the replace- 
ment 

without changing its value. Similarly, on the r.h.s. we can write 

d d 

F c (x 32 ) = -t—v c (x 23 ) = t—v c (x 32 ) (2.50) 
ox 2 ox 3 

The coordinate x 2 is now not involved in any operation, so for convenience it can be set 
equal to zero. Let us now proceed to transform the modified BGY equation into Fourier 
space. To do this we multiply both sides of the ^-independent form of (2.48a) by e tkxi 
and integrate over x\ with the condition |xi| > a. On the l.h.s. we have 

UZ + D dXietkXl (-^[ C '( a; i)- 2 ^V(a; 1 )]^, (2.51a) 
which after integration by parts becomes 

-2i sm(ka)C(a) + ik[C(k) - 2q 2 p] (2.516) 

On the r.h.s. we can extend the integration to the entire real line because this side 
vanishes for < a and remains finite at |xi| = a. Now it follows from (2.48b) that 
R(xi,x 2 ,x 3 ) is a fully truncated quantity and thus the integral involving R on the r.h.s. 
of (2.48a) is absolutely convergent. The integral over X\ of this term can therefore be 
done before the integral over x 3 . Doing this and using the convolution theorem in the 
first term gives that the Fourier transform of the r.h.s. equals 



-2Tpikv c {k)C{k) + Tq f°° dx 3 dVc } X ^ R(k; x 3 ) + qA(k) (2.52a) 

J— co UX 3 
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where, from (2.48b') and (2.48c) 



R(k;x 3 ) = £ s 3 p S2S M^)Q S2 sM^\k)-(j)e lkX3 C(k) (2.526) 



s 2 ,s 3 =± 

and 

qA(k) = 2iq 2 sm(ka)(p + -(a)+p- + (a))+qp- + (a) (<2+_(0, -cr\k) + Q-+(0, a\k)^ 

-qp + -(a) (g+_(0, a\0) + Q_ + (0, -a^)) (2.52c) 
(in deriving (2.52c) we have used the symmetry p + _(x) = p_ + (— x), where p SlS2 ( x ) '■ = 

To obtain the desired sum rule, we consider the leading order small-A; behaviour of 
(2.51) and (2.52). The term independent of k on both sides vanishes due to sum rules for 
the perfect screening of an internal charge: 

/oo 
dxC{x) = (2.53a) 
-oo 

an d _ 

Q SaSb (x a ,x b ;0) := / dxQ SaSb (x a , x b \x) = (2.536) 

J — oo 

which are expected to be true in all phases of Coulomb systems^ 17 ). The leading small-A; 
term is thus proportional to k. From (2.51), on the l.h.s. it is given by 

-2ikoC{o) - 2ikq 2 p (2.54a) 

while from (2.52) the leading small-A; terms on the r.h.s. are 

-2Tpikv c (k)C(k) + ikTq ^ s 3 [ dx 3 dv * Xs p S2S3 (0, x 3 )p S2S3 (0, x 3 ) 

s 2 ,s 3 =± J -°° ax 3 

+2iq 2 kalp + _(a) + P- + (a)} + ikqp-+(a)\p + -(0, -a) + p_ + (0, a)] 

-tkqp + _ (a) [p + _ (0, a) + p_ + (0, -a)} (2.546) 
where ^ 

dxxQ SlS2 (0,x 3 \x) (2.54c) 

-oo 

Let us now equate (2.54a) and (2.54b). The resulting equation can be simplified by 
recalling that p ++ (x) — for all \x\ < 2a, and so 

C(*) = -q 2 (p + _((7)+p_ + ((7)) 

Furthermore, v c (k) = n/\k\, and by symmetry p ++ (x) = p (x) = for all x in any 

phase. We thus obtain the sum rule 



lim7r^C'(fc)/|A;| 



fc->0 



1+ 2^/ dx ^(P-+( x "> p -+( x "> ~ P+-( X )P+-( X )) 



2p 
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+ J~P (P-+W p -+W ~ P+-(°)P+-(°)) ( 2 - 55 ) 

which is to be obeyed in all phases. However we have commented in subsection 1.1 that 
for the charge ordered system C(r) is expected to exhibit a decay which is o(l/r 2 ), except 
possibly on the phase boundary. Such a decay implies C(k) = o(\k\) and thus the l.h.s. 
of (2.55) vanishes, leaving us with a simplified form of the sum rule which is to be obeyed 
by the system in all phases except possibly on the phase boundary. Note in particular 
that the dipole moment p + ~(x) is non-zero in both phases. In Ref. [9] it was claimed that 
p + -(x) as calculated from the exact results at T = 1 vanished. However, in reviewing that 
calculation, an error has been found and the correct conclusion is that p + _(x) is non-zero 
at T = 1, in agreement with the general property. 

2.4.1 Verification of the sum rule for the T = reference system 

We have shown in subsection 2.1 that the T = reference system perfectly screens an 
internal charge. Since the perfect screening properties (2.53) are the two fundamental 
assumptions in deriving the sum rule (2.55) from the BGY equation (2.48a), we therefore 
expect the sum rule to be satisfied in the reference system. Let us check this property. 

The reference system of subsection 2.1 consists only of point particles, so the hard 
core width in (2.55) needs to be taken to zero. Also, the final terms on the r.h.s can 
be simplified due to the symmetries p + _(x) = p_ + (x), p + _(x) = —p_ + (x), and the pair 
potential in the reference system is identically zero. In these circumstances the sum rule 
(2.55) reads 

= 1-V(0K_(0 + ) (2.56) 

(p H (x) is discontinuous at the origin so it is necessary to specify its value on one side). 

From (2.9), p H (0) = 2p 2 . It remains to calculate p ^(x). From the definitions (2.54c) 

and (2.46c) we must first calculate p+L + (xi, x 2 , x 3 ) and p+!__(xi, x 2 , x 3 ), which can be 
done by an appropriate extension of the method given in subsection 2.1 to calculate 
p ( +l(x). We find 



( \ = l P 3 ^ 1 + e- 4 P\ X3 ~ X2 \ + e- 4 ^ 2 -* 1 ! + e -M^-x!\^ for Xl < ^ < Xs 

P+-+[xi, x 2 , x 3 ) - | ^ + ^i^-^i _ e -4 P \ X2 - Xl \ _ g-^s-sil), for Xl<X3< X2 



(2.57) 



and 



( ( p3(l- e -4p\*3-*2\ +e -Mx2-*i\- e -Mx3-*i\), forx 1 <x 2 <x 3 

P+—{X 1 ,X 2 ,X 3 ) - { p 3( 1 + e -4 p \x 3 -x 2 \ +e -4p\x 2 -x 1 \ +e -4p\x 3 -x 1 \^ for ^ < ^ < ^ 

(2.58) 

Using (2.57), (2.58) and (2.9) to form Q + _(xi, x 2 \x) as given by (2.46c) we can check 
the perfect screening property 

/oo 
dxQ + _(x 1: x 2 \x) = (2.59) 
-oo 

which together with (2.53a) is a fundamental assumption in the derivation of (2.53). We 
can also calculate the first moment of to obtain 

P+ _(^K-W = ±fe- 4 ^ (2.60) 
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where the positive sign is to be taken for x > 0, while the negative sign holds for x < 0. 

From (2.60) we see that indeed the sum rule (2.55) is obeyed for the T = reference 
system. 

3. THE SYSTEM WITHOUT CHARGE ORDERING 
3.1 Decay of correlations in the conductive phase 

For T < 2 the two-component log-gas without charge ordering is expected to perfectly 
screen an infinitesimal charge density. As noted in subsection 1.2, a linear response 
argument*- 8 -* then leads to the sum rule (1.4b) for the asymptotic behaviour of the charge- 
charge correlation. In this subsection we will note another derivation of (1.4b), using a new 
method due to Jancovici^ 18 ^ , which is based on a linear response argument and macroscopic 
electrostatics. In ref. [18] the conductive phase is characterised by the applicability of 
macroscopic electrostatics at large length scales in the system. The log-gas in a conductive 
phase is thus considered as an infinite metal line obeying the laws of two-dimensional 
electrostatics. When combined with a linear response relation, describing the change 
in the potential due to the addition of an external charge, this characterisation implies 
(1.4b). Furthermore it is clear that (1.4b) only applies to the leading non-oscillatory term 
of the charge-charge correlation only: the use of macroscopic electrostatics implies the 
charge density must be smoothed over some microscopic distance and thus the oscillatory 
terms averaged to zero. 

3.1.1 Abe-Meeron type resummations 

In three-dimensional classical Coulomb systems a systematic way to study the high- 
temperature conductive phase is via the Abe-Meeron^ 12 ) diagrammatic expansion of the 
Ursell function h SlS2 (r) (:= pJ lS2 (r) / p 2 ) . This approach is also applicable to the present 
system. 

Briefly, let us recall (see ref. [19] for a detailed recent review of the method, and an 
extension to quantum systems) that h SlS2 (r) is given by the sum of all the Mayer graphs 
built with the Mayer bonds 

f SlS2 (r)=e- Tv ^(r)-l. 

(here the potential v SlS2 corresponds to the regularized logarithmic Coulomb potential 
between unit charges of signs si and s^). Once the chain summations have been performed, 
the above set of Mayer graphs is exactly transformed into an a new set of 'prototype' 
graphs II with the same topological structure and two kinds of resummed bonds: 

(i) the Debye-like bond SlS2 (r) defined as the sum of all the convolution chains built with 
the Coulomb potential siS2V c (r), and such that 

1 (us _ T Sl s 2 v c (k) . 

(convolutions of this bond are to be excluded). 

(ii) the resummed bond /#: 

f R (r) = e - r («'i. 2 W-i- 2 «c(r))+*. 1 . 2(r) _ 1 _ 
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At large-distances, v SlS2 (r) — sis 2 v c (r) vanishes and so /^(r) behaves as 2 iS2 (r)/2. 
Now, since v c (k) ~ n/\k\ for small \k\ we deduce from the definition of the Debye-like 
bond that 4> SlS2 { r ) decays as 1/r 2 for large-r. Therefore all the resummed bonds (i) decay 
at least as 1/r 2 while the resummed bonds (ii) decay as 1/r 4 ; the Ursell functions h slS2 (r) 
therefore exhibit a 1/r 2 decay in the conductive phase. This is to be contrasted to the 1/r 4 
decay found in the high temperature phase of the system with charge ordering, where the 
Abe-Meeron diagrammatics do not apply due to the presence of a many body potential 
inducing the charge ordering constraint. 

3.2 Decay of correlations in the dipole phase 

Our objective in this subsection is to use low fugacity expansions to deduce the large- 
r decay of the truncated two-particle distributions and the charge-charge correlation for 

T > 2. Since in the system without charge ordering p ++ (r) = p (r) and p + _(r) = p~+(r) 

it suffices to consider p+ + (r) and p+_(r). We will study in detail the terms of 0(C 4 ) since, 
as is argued below, the same asymptotic decay is expected of terms of higher order in (. 

At order ( A the low fugacity expansions of the two particle correlations can be cal- 
culated from (2.21), (2.22) and the explicit forms of the partition functions which occur 
therein. We find 



p+ { +{r) = ( 4 J dx x dx 2 



r(x 2 - x-i) 



(x 2 — r)(xi — r)x\X 2 



1 



bi| r |r 



x 2 \ 



(3.1) 



and 



p?-V) = 

dx\dx 2 



xi(x 2 - r) 



r r bi 



x 2 



where 



r(x 2 - x 1 )(x 1 - r)x 2 
1 

S/(zi, x 2 ) := - [f(x 1 ,x 2 ) + f(x 2 , x\, 



-S- 



and the notation * denotes the region of integration 

{(xi, x 2 ) G R 2 : \xi — x 2 \ > a, \x±\, \x 2 \ > a, \x 1 

From these expansions, and (2.24), we have 
C< 4 »(r) = 



\x 2 — r\ > a} 



(3.2) 

(3.3a) 
(3.36) 



dx\dx 2 



r(x 2 - xi) 



(x 2 - r)(xi - r)x x x 2 



2S 



xi(x 2 - r) 



r(x 2 -x 1 )(x 1 -r)x 2 



+ 



r r \xi — x 2 \ T 



(3-4) 

In Appendix B a method of analysis of the large-r behaviour of the integral (3.4) will 
be given. We find that for 2 < V < 4 the region of the integrand which gives the leading 
order contribution comes from the physical configurations a of Figure 5 in Appendix B, 
while for T > 4 the configurations (3 of Figure 5 in Appendix B give the leading order 
behaviour. Furthermore this remains true of the integrals (3.1) and (3.2). Here we will 
calculate the contribution to the large-r behaviour of the integrals in (3.1) and (3.2) from 
these configurations. 
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For the configurations a of Fig. 5, which give the leading large-r behaviour for 2 < 
T < 4, we can replace |xi — r| r and |x 2 — r| r by r r in (3.1) and (3.2), provided we 
also multiply by a factor of 2 due to the invariance of the integrand under the mappings 
x 2 — r \— > x 2 , Xi — r i— > x±. This gives 



p+y{r) 



2CV 2 



J dx\dx 2 



1 


x 2 




Xi 


r 


2 


Xi 


r 


x 2 


r 



Fi 



and 



P?- 4) (0 



2CV 2 



(ar) 1 







r 






1 




1^2 - 


X 1 


r 


^2 


r 


Xi - x 2 


r 



dx\dx 2 

where * now denotes the region of integration 

{(xi,x 2 ) G R : |xi - x 2 | > 1, |xi], |x 2 | > 1} 



Xi 



(3.5) 



(3.6) 



(3.7) 



and we have removed the a dependence from the integrals by changing variables X\ \— > 
o"Xi,x 2 I— > <7X 2 . Note that the integrals in (3.5) and (3.6) are conditionally convergent: 
for fixed X\ (or x 2 ) one must add the contributions of positive and negative x' 2 s (or x[s). 

For T > 4, when the configurations (3 of Fig. 5 give the dominant contribution to the 
asymptotics, the integrals (3.1) and (3.2) are analysed by expanding the integrands for 
Xi near and x 2 near r (again we need to multiply by a factor of 2 to account for the 
symmetry in the interchange of x\ and x 2 ). This gives 



and 

where 

A(u,v;r): , 



P+-(r) 



* dudv 



\uv\ 



dudv 



\uv\ 



TA(u, v; r) + 



r 2 fuv\ 2 



—TA(u, v; r) + 



r 2 /uv\ 2 

T \7 2 



(3.8) 



(3.9) 



uv 1 
, ■ + 3 



u — V 



+ 





3 




3" 


1 






I) 




(f) 




4 





-■9- 



and * denotes the integration region 

{(u,v) G R 2 : \u\ > a, \v\ > a} 



3 

(3. 



Oa) 



(3.106) 



In performing the integrations, only those terms in the integrand even in u and v survive, 
so (3.8) and (3.9) simplify to give 



and 



1 



(r-3)a r - 



(r - 3)a r - 3 



(3.11) 



(3.12) 



The formulas (3.8) and (3.9) have some similarities and some differences in relation to 
the analogous results for the 2dCG (ref. 7 eq. (3.17)). One similarity is that the 0(l/r 4 ) 
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terms proportional to T 2 result from squaring the dipole-dipole potential uv/r 2 between 
the neutral clusters which form the dominant configurations. Another similarity is that 
A(u, v;r) is the multipole expansion (appropriately truncated) of the potential between 
the neutral clusters, and that all terms odd in the mappings v i— > — v or u i— > —u vanish. 
A crucial difference is that in one space dimension the log-potential is not harmonic, so 
unlike the two-dimensional case, the integration over A(u, v; r) does not vanish. The first 
non-vanishing term is the quadrupole-quadrupole potential, which has the same 0(l/r 4 ) 
decay as the square of the dipole-dipole interaction. This term gives a contribution of the 
same magnitude to both (3.8) and (3.9) but of opposite sign. 

At the coupling T = 4 the leading order behaviour of p^.+\r) is given by the sum of 
(3.5) and (3.11), and the leading order behaviour of p+^(r) is given by the sum of (3.6) 
and (3.12). At the coupling T = 2 all regions of integration analysed in Appendix B were 
bounded by an 0(l/r 2 ) decay. However, it is possible to evaluate the integrals in (3.1) 
and (3.2) directly by using partial fractions. This is done below, where it is shown 

AV) ~ ^ (3.13) 

and 

~ (3.14) 

In fact these results are precisely the values given by (3.5) and (3.6) at T = 2; the 
coefficient of the l/r r term is therefore continuous as T — > 2 + . 

For general n > 2 the large-r behaviour of p^ n \r) and p^L n \r) can be deduced 
by hypothesizing the generalization of the configurations a and (3 of Fig. 5, which give 
the leading large-r contribution, for 2 < T < 4 and T > 4 respectively, to the integral 
expression for p++\r) and p+^(r). The obvious generalization of configurations a, which 
should give the leading large-r expansion for 2 < T < 4, have n — 2 — p mobile particles 
about one root particle and the remaining n + p, with < p < n — 2 mobile particles 
about the other root particle. The mobile particles are to be distributed so that the total 
charge of the cluster about and including one root particle is +q while the total charge 
about and including the other root particle is —q. At large distance the effective potential 
is thus q 2 log r which gives the required behaviour as 

P T\r) - (3.15) 

pT\r) ~ ^^C 2 " (3.16) 

This is the same order decay as exhibited by (3.5) and (3.6) in the case n — 2, and 
furthermore the amplitudes could, in principal, be given integral representations analogous 
to those for n = 2. 

In particular, the coefficients in (3.15) and (3.16) should be finite as T — > 2 + . This 
would imply that C(r) as calculated from a resummation of (3.15) and (3.16) would exhibit 
the same power law decay 0(l/r r ), independent of the fugacity (. We therefore expect 
that the critical line separating the conductive and insulator phases to be at T = 2, 
independent of the fugacity, which is contrary to what we have found for the charge 



and 
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ordered system, but in agreement with previous studies^ 3 ' 20 ) . Thus the nested dipoles of 
subsection 2.2.3 do not give the dominant contribution to C(r) in the system without 
charge ordering. In fact the analysis in Appendix B shows that the leading contribution 
of nested dipoles cancels out when the charge ordering constraint is removed. 

The generalizations of the configurations (3 of Fig. 5, which should give the leading 
behaviour for T > 4, are neutral clusters about and including the root particles. The 
potential V 2 n between these neutral clusters can be expanded as 

V 2N = Wo + W r + U 0r 

where W and W r are the electrostatic energies of the neutral clusters about and r 
respectively. For large-r, the potential Uq t can be written as a multipole expansion in r. 
Analogous to the case n = 1, there will be two classes of terms contributing to the final 
leading order expansions 

pVi"\r) ~ C A :; {T) (3.17) 



r 

and 

T(2n) 

, . . 

r 



^» ~ C2Mg 4 - (r) (3.18) 



3.2.1 Decay at Y = 2 

Here we will show how (3.13) and (3.14) can be derived. Consider for definiteness (3.13). 
From the Cauchy identity 



r(x 2 - xi) 



det 



y (x 2 — r)(xi — r)xix 2 J 
we see that at T = 2 (3.1) can be simplified to read 

p++\r) = -C 4 J dxidx 2 



l/x\ l/(xi — r) 
l/x 2 \/{x 2 -r) 



1 



(x 2 - r)(xi - r)x\X 2 



(3.19) 



where the integration domain * is that specified by (3.3b). The integration over x 2 can 
now be performed using the simple partial fraction expansion 



(x 2 - r)x 2 r 



— --) 

x 2 -r x 2 J 



(3.20) 



We must consider separately the regions — oo < x\ < —2a, —2a < x\ < a, a < x\ < 2a 
and 2a < Xi < 1/2 (for Xi > 1/2 an identical contribution to (3.19) is obtained). This 
gives, after changing variables x x \— > ru, 



2( 4 



/— 2a/r r 
+ / 



V2 du 



2a/r [U — 1)U 



log 



u + a/r 



u — a/r 



+ 



+ 



- a / r du 



2a/r (u — l)u 
2<y I v du 
air (u— l)u 



log 


a/r 


- log 


-1 + a/r 


u — a/r 


u — 1 — a/r 







log 



u + a/r 



a/r 



-log 



u — 1 + a/r 



— 1 — a/r 



(3.21) 
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where we have ignored an additive term log[(l + o jr)j(\ — ojr)\ in each integrand since 
it decays as r — > oo. In fact the only portion of the integrand which does not decay as 
r — > oo comes form the region u = 0(a/r). Setting u = xa/r and keeping only those 
terms which are non-zero gives 

p T + f (r) ~ ^ (/ 2 dx I log(l + x) + cfa 1 log |±i) (3.22) 
Evaluating these integrals gives the result (3.13). 
3.3 A sum rule from the BGY equation for C(x) 

Our objective in this subsection is to derive sum rules analogous to (2.55) for the system 
without charge ordering. We will consider first the case when the logarithmic attraction 
between opposite charges is smoothly regularized. In this case we will provide an illus- 
tration of the sum rule using a solvable "parallel line" model at V — 2. The other case to 
be considered is the hard core regularization used in the above subsection. 

Let us denote the smoothly regularized potentials by v SlS2 (x) and define the corre- 
sponding force by 

F sm(x) = ~^ Vs ^( x ) ( 3 - 23 ) 

(using this notation rather than that of the subsection 2.4 allows for the possibility of 
different regularization between like and unlike species). The BGY equation for p SlS2 ( x i2) 
is then given by (2.46a) with SiSjF(xij) therein replaced by F s . s .(xij). Assuming the 
symmetry 

F SlS2 (x) = F ( _ Sl)( _ S2 )(x), 

it is straightforward to show from the BGY equation for p SlS2 (xi 2 ) that the BGY equation 
for the charge-charge correlation in this system is (recall the derivation of (2.48a) from 
(2.46a)) 

° [C(x 12 )-2q 2 pS(x 12 )} 



dx 



/oo 
dx 3 [F ++ (x 32 ) - F + ^x 32 )] C(x 13 ) 
-oo 

/oo 
dx 3 F ++ (x 32 ) lp++(x 32 )Q ++ (x 2 , x 3 \xi) - (p/2q)C(x 31 )} 
-oo 
roo 

-2Tq dx.F+.ix^ip+^x^Q^x^x^-ip^Cix^)] (3.25) 



We now follow the procedure detailed in subsection 2.4 to study the small-/c behaviour 
of the Fourier transformed version of (3.24). The final result is the sum rule 

lim Tt (3C {k) / \k\ = 1 - — f° dx dv+ ~^ p„ + (x)p„ + (x) (3.25) 

k— »0 p J— oo OX 

where the dipole moment p_ + (x) is defined by (2.54c). In the distinct phases different 
terms in (3.25) have special values which allows further simplification. In the conductor 
phase T < 2 we expect p_ + (x) = 0, which implies 

ir/3C(k)/\k\ = 1 (3.26) 
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This result is equivalent to (1.4b). On the other hand, in the insulator phase (r > 2) the 
expected decay (1.4a) implies the l.h.s. of (3.25) vanishes and thus the dipole moment 
P-+(x) is non-zero. 

In the case of a hard core regularisation of the logarithmic potential, working analogous 
to that given in subsection 2.4 for the charge ordered system shows that the sum rule (2.55) 
is still applicable. Without charge ordering the correlations obey 



P+-( X )P+-( X ) = ~P-+( X )P-+( X ) 



so the sum rule reads 



lim ndC(k)/\k\ — 1 — — / dx ^-p + ^(x)p + 

k^O p 7_oo dx 



qp 



-P+-( <r )P+-( tr ) 



(3.27) 



In the insulator phase (r > 2) we expect the behaviour (1.4a) so the l.h.s. of (3.27) will 
vanish. In the conductor phase (r < 2) the dipole moment p+~(x) should vanish and 
thus (3.26) is reclaimed. 

3.3.1 Verification of the sum rule for a parallel line model at T = 2 

The 2dCG is exactly solvable at the coupling T = 2^ u \ where the multi-particle distri- 
bution functions can be calculated. The method of exact calculation used in ref.[21] can 
be adapted to a model in which the negative charges (coordinates Xj,j = 1, . . . , N) are 
confined to a line and the positive charges (coordinates Xj, j = 1, . . . , N) are confined to 
another line, parallel to the first line and separated by a distance e. The pair potential 
between opposite charges of unit strength is thus given by 



v _ + ( x - X) = log[(x - Xf + e 2 ] 1/2 



(3.28) 



while the pair potential between like charges is the Coulomb logarithmic potential. Notice 
that (3.28) is a smoothly regularized Coulomb potential. As such, the system should obey 
the sum rule (3.25). Let us verify this fact at the solvable coupling T = 2. 

To verify (3.25) we need to evaluate p- + (xi,X2), which requires the two and three 
particle distributions. Using the technique of ref.[21], we can easily evaluate the general 
distribution function as 



P+...H — ...-(Ai, . . . X m ; Xi, . . . , x n2 ) 



where 



and 



det 



[G ++ (Xj, -^fe)]j,fe=l,...,n 

[G+-(Xj, Xk)]j=l,-,ni 
k=l,...,n2 



G ++ (x) = G (x) 



G_Ax) = G+Jx) 



[G-+(Xj, X k )\ 2=l,...,n 2 
[G (Xj, Xk)]j,k=l,...,n 2 



= /„ d ''T 



+ (l/27rC) 2 e 4 ^ 7 



o?7 



2tt< 1 + (l/2vrC) 2 e 4 ^ 



(3.29a) 
(3.29b) 

(3.29c) 
(3.29d) 
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We want to use these formulas to evaluate p + ~(x) as given by (2.54c). For this purpose 
we note from (2.46c) that 

P- + (xi,x 2 )Q-+(xi,x 2 ; x) = q (p+_ + (x, x 1 ,x 2 ) - p T _ + (x, x 1 ,x 2 ) + a(x; xi,x 2 fj (3.30a) 
where 

a(x;x 1 ,x 2 ) = pp+_(x,xi) + pp+ + (x,x 2 ) - pp T _(x,x 1 ) - pp T + (x,x 2 ) 

+ (5(x - x 2 ) - S(x - x±)) p- + (xi,x 2 ) (3.306) 

Now it follows from the perfect screening sum rule (2.53a) for a single internal charge 
(which can be verified using (3.29)), and the dependence of the truncated two particle 
distributions pJ 1S2 (x,x') on \x — x'\, that 

/oo 
dxxa(x;xi,x 2 ) = q(x 2 - xi)p^ + (xi,x 2 ) (3.31) 
-oo 

Next, to evaluate 

/OO y . 

dxx (p+_ + (x,x 1 ,x 2 ) - p^_ + {x,xi,x 2 )\ (3.32) 
-oo v ' 

we note from (3.29) that 

p T + _ + (x,x u x 2 ) - pT_ + (x,X!,x 2 ) 

= 2Re (G ++ (x, x 2 )G+-(x 2 , x{)G- + (x u x) - G ++ (x,x 1 )G + -{x 1 ,x 2 )G- + {x 2 ,x)) 

(3.33) 

which shows it suffices to evaluate 

/oo 
dxxG ++ (r, x 2 )G- + (xi, x) (3.34) 
-oo 

This latter task can be carried out with the aid of the formula 

/oo 1 f) 
dxxe 2mx ^-^ = ^— 5(7i - 72) (3-35) 
-00 2m c>7i 

and integration by parts. Once having evaluated (3.32), we substitute the result in the 
integral over r of (3.30a) together with (3.31) to conclude 

p- + (xi,x 2 )p-+(xi,x 2 ) = 

coo p 2-Ke^+2m^(x 2 -x 1 ) roo 2ire^ 1 +2m^ 1 (x2-x 1 ) 



iqe 



(27TC) 2 



Jo dl 1 + (l/27rC) 2 e 4 -7 I d ^ (TTW^CP^F ~ " X2 \ 



+q(x 2 - x 1 )pT + (x 1 ,x 2 ) (3.36) 
With the result (3.36) and the formula 

dv_ + (x) 1 ( 1 1 



dx 2 \x + ie x — ie 
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for the force between opposite charges, all quantities in the integrand on the r.h.s. of the 
sum rule (3.25) are known. The integral can now be carried out by aid of the formula 

3 27ri(7i— 72)2; 



_ 0~„-„27r(7i-7 2 )e w J > '1 '2 

0, otherwise 



/>p _ +Wp _ +w ^M. f I-*™ (3,7) 



/■oo „z7rn7i-7 2 ;x ,• 

/ rfx ; = -2nie 27T ^ ,2)€ x 

J —oo X ~\~ 16 I 

We find 

r ^ n .M„ .U ^+ix) = qp(^_ 

1 + (1/27TC) 

which thus evaluates the r.h.s of (3.25). 

To evaluate the l.h.s. of (3.25), we note from the exact formula (3.29) the large-x 
expansion 

r( x _ 1 

[X) ~ 2{7ixy 1 + 1/(2ttC) 2 

which implies 

Substituting (3.37) and (3.38) in (3.25) we see that the sum rule is indeed valid. 

3.4 Decay of 3 and 4 body correlations in the dipole phase 

The decay of the higher order correlations can be studied using the low fugacity expansion 
method of subsection 3.2. Alternatively, the decay of the 3 body correlations can be 
deduced from knowledge of the decay of the 2 body correlations and use of the BGY 
equation. Let us consider this latter method. 

Suppose for defmiteness that the logarithmic potential is smoothly regularized at the 
origin and v slS2 = SiS 2 v, so that the BGY equation (2.46a) applies. Making the changes 
(2.49) and (2.50) allows the Fourier transform with respect to x\ of both sides of (2.46a) 
to be taken, with the result 



fo(x< 
dx 3 



+ s 2 rj_ 



0,x 3 ) -# iaa _(M,s 3 )) ( 3 - 39 ) 



Now from (3.5) and (3.6), for 2 < T < 4 

PlsM ^ j~jr ^d so plJk) ^ \k\^ (3.40) 

(here and below the numerical amplitudes in the asymptotic formulas are omitted). The 
first term on both sides of (3.39) is thus 0(|/c| r ) (here we have also used F(x) ~ 1/x). 
Consider the second term on the r.h.s.. Since v c (k) = n/\k\, non-analytic behaviour 
results from both the leading non-analytic and leading analytic term in the expansion of 
C{k): 

C{k) ~ Ikf- 1 + k 2 (3.41) 

This term is of a lower order in k than the previous two considered and so must be 
balanced by the 3 body term: 



° dx 3 ^^(pl S2+ (k,0,x 3 )-pl S2 _(k,0,x 3 )) ~^ pi nsgn(k)(\k\^ + k 2 ) (3.42) 



dx 
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Note that for T > 3 the 0(k 2 ) term is dominant. Taking the inverse transform suggests 
that for large- 1 xi | and fixed x 3 

!q / sgn(xi) \ 2 < r < 3 
o (vf r > 3 (3,43) 

At r = 3, an extra logarithmic factor is expected. 

The non-perturbative result (3.43) is in agreement with a low fugacity analysis. In 
particular, the 0{1/ x\) behaviour of pj_ + (xi, 0, x 3 ) at 0(( 4 ), which is dominant for T > 3, 
originates from the configuration in which the mobile particle of charge — s forms a dipole 
with the root charge at x\. Explicitly, the first term of the multipole expansion of this 
charge configuration which gives a non-zero contribution to the cluster integral is 

d 2 d 1 
(u - x 1 ) 2 (x 3 - x 2)-^i-^r (- lo £ - ^1) ~ ^3 (3.44) 

For the 4 body correlations, an analysis of the BGY equation for the 3 body correla- 
tions analogous to that given above leads to the conclusion that for all T > 2 

P«(-, 1 )«(-«)( a; i' x 2, x 3 , x A ) ~ ^ (3.45) 
as \xi — x 3 \ — > oo with \xi — x 2 \ and \x 3 — x±\ fixed. 
3.4.1 Exact decay of a three particle correlation at Y = 4 

A version of the two-component log-gas system without charge ordering, in which the 
two species of charges are confined to separate, interlacing sublattices (lattice spacing r), 
is exactly solvable at the isotherm T = A^ 10 \ In particular, the exact expression for the 
three particle correlation involving opposite charges is ^ 

^__(n,mi-l/2,m 2 -l/2) 

= ^Tr [B+_(n,mi - l/2)B__( mi - 1/2, m 2 - l/2)5_ + (m 2 - 1/2, n)} 

+^Tr [B + _(n, m 2 - l/2)B__(m 2 - 1/2, m x - l/2)B_ + {m 1 - 1/2, n)\ (3.46) 

where 

B + .(»,m. - V2) = ^ M) , B_ + (m 2 - 1/2, „) = ^ 



with 



wo =—(i+o r 

r Jo 



p2irilt 



D(t) 



WO = -^(1-0^*^(1/2-0 



Mi ) = ^f(l + «/ o *^(t-l) 
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D(t) :=l-lQt(t 2 -t+l/8)+£ 2 , £:=C 2 7r 4 /T 4 , Zi := n-(mi-l/2) / 2 :=m 2 -l/2-n 
and 

B__ (mi - 1/2,^-1/2) = (M« 

with 

/M0 = / cft^yK 2 -^(l-t)-^-l/2) 2 ] 

„1 2-wilt 

mi) = ^AKi y o dt_(t- 1/2) 

/MO = ^/>gJ(«-i/W-i) 

/ 3 := mi - m 2 

We seek the asymptotic behaviour of (3.46) for |ra2| — > oo with n and mi fixed. Now, 
straightforward integration by parts of the above integrals gives 

As(z)~7 /3 3 i(z) ~ ^ /M0~^ 
/M0~j^ /M0~j Ai(0~j^ 



and thus 

B__(mi - l/2,77i2 - l/2) J B_ + (m 2 - 1/2, n) ~ I «J5 | (3.47a) 



i i 

m2 m5 



B + _(n,m 2 - l/2)B__(m 2 - 1/2, mi - 1/2) ~ [ 7 7 | (3.476) 
Substituting (3.47) in (3.46) gives that 

p^__(n,mi - l/2,m2-l/2) 
in agreement with the prediction (3.43). 



m 2 ^oo m 3 



4. SUMMARY AND COMMENTS 

We have undertaken a systematic study of properties of the particle and charge correlation 
functions in the two-dimensional Coulomb gas confined to a one-dimensional domain, with 
and without charge ordering. As a result some new properties have emerged. 

Consider first the system with charge ordering. In the high temperature phase we have 
found an 0(l/x 4 ) decay of the two-particle correlations as given by (2.18). The system 
does not screen an external charge density, as this would lead to (1.4b) and thus an 
0(l/x 2 ) decay. However the underlying reference system at T = does screen an internal 
charge. This allows an expansion of the large- a; two-particle correlations to made about 
r = and leads to the 0{l/x 4 ) decay. For T -> 2+, it was shown in subsection 2.2 that 
the asymptotic charge-charge correlation can be resummed to all orders in the fugacity, 
and that the Kosterlitz renormalization equations result. The radius of convergence of 
the quantity A (2.19) as given by (2.37) gives the dependence of the critical coupling 
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on the fugacity, analogous to the situation with the resummed dielectric constant in 
the 2dCG^ 7 \ Furthermore, the configurations contributing to the asymptotic charge- 
charge correlation were identified as nested dipoles, analogous to situation in the scaling 
region of the Kosterlitz-Thouless transition of the 2dCG. Away from criticality in the low 
temperature phase the two-particle correlations have an 0(l/r 2 ) decay as given by (2.44). 
In a low fugacity perturbation expansion this behaviour can be understood as resulting 
from the dipole-dipole interaction between neutral clusters of particles about each root 
particle. In all phases the correlations must obey the sum rule (2.55) involving the dipole 
moment p + ^(x) of the charge distribution induced by two fixed opposite charges. This 
sum rule is analogous to the sum rule (1.3b) for the 2dCG, which relates the dielectric 
constant to p- + . 

Now consider the system without charge ordering. In the high temperature phase 
T < 2 the system is conductive and the charge correlation must exhibit the asymptotic 
behaviour (1.4b). The low fugacity analysis of the two-particle correlations given in 
subsection 3.2 imply a 0(l/x r ) decay for 2 < T < 4 and an 0(l/x 4 ) decay for T > 4. 
In the context of the low fugacity expansion this behaviour results from configurations 
consisting of clusters about the root charges which have a charge imbalance equivalent to 
that of one particle dominating for 2 < T < 4, while neutral clusters dominate for T > 4. 
Furthermore, for r — > 2 + the coefficients of the low fugacity expansion of the large- X two- 
particle correlations are all finite; there is no dominance of nested dipole configurations 
which in fact cancel without charge ordering. The finiteness of the coefficients indicates 
that the phase transition occurs at T = 2 independent of the fugacity. With a smoothly 
regularized potential, we have presented the sum rule (3.25) involving the dipole moment 
p H (x), and explicitly verified it on a solvable model at the coupling T = 2. 

As far as large-distance behaviours of the internal correlations are concerned, the 
system with charge ordering is very similar to the familiar 2D Coulomb gas. Indeed, 
the charge correlations have a "fast" decay in the high temperature phase (l/|a;| 4 in 
comparison to an exponential decay in 2D), and as a density-dependent power law in the 
low-temperature phase (l/|x| rA in comparison to l/|r| r / £ in 2D). The partial screening of a 
given pair (with sizes X in ID and R in 2D) by smaller pairs takes the same mathematical 
form for both systems, because of the occurence of the integrals 

f d 

/ dx — log Id (4.1a) 

Ja<x<x dx 

and 

/ dffviogr) (4.16) 

Ja<r<R V ' 

which both diverge logarithmically when X and R go to infinity. The integral (4.1a) 
arises from the dipole-dipole potential between the screening pair and the root charges. 
The resulting screening contribution is proportional to the average of the dipole carried 
by this pair, which does not vanish by virtue of the charge ordering constraint. 

Of course, the screening contribution associated with (4.1a) disappears in the system 
without charge ordering. The leading contribution of a screening pair then involves 



/. 



dx (^-log\x\\ (4.2) 

a<x<X \dx J 

which is linked to the fluctuations of the dipole-charge potential between this pair and the 
root charges. The integral (4.2) does converge in the limit X — > oo, and consequently the 
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collective effects do not affect the power which controls the algebraic decay of the charge 
correlations. The Kosterlitz-Thouless transition temperature is then expected to remain 
constant (T = 2) at low densities. This is contrary to the findings of a recent computer 
simulation^ 23 ), where the transition coupling appeared to decrease with decreasing density. 
However, we believe these findings are linked to finite-effects which become more and 
more important as the density decreases. Indeed, in the zero-density limit and with 
1 < T < 2, the proportion of free charges goes to zero (see e.g. ref. [15]), a behaviour 
which is reminiscent of the short-range collapse at T = 1 for the system without hard-core 
regularization. In this regime of parameters, the finite systems considered in the numerical 
simulations are not efficient for perfectly screening external charges because of the very 
small number of free charges which are available. A similar mechanism also occurs in 
2D( 24 \ and even in 3D^ 25 ) for the restricted primitive model with 1/r interactions (in that 
case there is no Kosterlitz-Thouless transition). Eventually, we note that at high densities, 
despite difficult problems of sampling, the numerical simulations'- 23 -' are compatible with 
the theoretical predictions given here and in ref. [3]. 

Our predictions should be valid at low densities. At higher densities, the critical 
Kosterlitz-Thouless transition line in the (p, T) plane might bifurcate, at a tricritical 
point, into a first-order liquid-gas coexistence curve. In fact, for the 2D Coulomb gas, 
this has been observed in computer simulations by Caillol and Leveque*- 26 -* . More recently, 
an extended Debye-Huckel theory of Levin et al/ 27 -* suggests a tricritical point at a much 
lower density. To our knowledge, such first-order transitions have not yet been observed 
for the present models in ID. At a theoretical level, their study is beyond the scope and 
the methods of this paper. In particular, a correct description of eventual oscillatory 
correlations^ 28 -* beyond a Fisher- Widom line requires further resummations of the low 
fugacity expansions combined to a detailed account of the short-range effects which depend 
on the hard-core regularization. At the moment, note that for the model without charge 
ordering, the exact results for the solvable isotherm T = 2( 10,11 ) seem to exclude the above 
bifurcation process. 
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Appendix A 



Our objective in this appendix is to show how the integrals in (2.23) can be analysed 
to determine the asymptotic charge density C^\r). We recall C^\r) is defined as the 
portion of C(r) proportional to ( 4 which contributes to the correct leading order behaviour 
of (2.19) in the limit T -> 2+. 

We will illustrate our method by calculating the contribution to C^\r) from p+ + (r). 
Scaling the integration variables by 

yi i-> ry x and y 2 i-> ry 2 

we read off from (2.23a) that 



P++( r ) - £-2 



•1-a/r poo ( y 2 — y 1 \ r 2r-2 



rL—a/r roo 

/ dyi / dy 2 

Ja/r Jl+a/r 



a/r Jl+a/r \ (y 2 ~ l)y 2 (l ~ 2/l)j/l / (r - l) 5 



(Al) 



The fixed root points are now at and 1. Let /i be a positive constant such that a/r < 
|i « 1. Then for the region of integration p < y x < 1 — p and y 2 > 1 + // in (Al) 
the integrand is bounded uniformly in r and is integrable. Hence the contribution to 
the double integral is 0(1), and from (2.19) the corresponding contribution to A^ 4 ) is 
0(£ 4 /(r — 2)). To analyse the contribution from the remaining region of integration we 
decompose it into the subregions indicated by the particle configurations in Fig. 2. 
Consider subregion A in which 

a/r < y± < /i and a/r < y 2 — 1 < /i 

Since y\ and y 2 — 1 are small variables we can expand the integrand: 

f 7 TTTT S ) ~ 7 + r V^ + ■) 

V (2/2 - 1)2/2(1 - 1/1)2/1 / (2/2 - l) r 2/i 

This shows that the leading large-r contribution to the double integral from subregion A 
is 

W^r + (f^" 2 " 1)2 (A2) 

Substituting (Al) in (A2) shows that the first term cancels. To leading order for large-r 
the second term gives a contribution 

t(4)/ x C r 

P++ (r) ' 



(r - 2)V 2 



The first moment of this term (considered as a function defined on [a, oo))is infinite. 
However the entire second term in (A2) is exactly cancelled in the formula (2.24) for C{r) 
by an identical term contributing to p+l T (r), so there is no net contribution to . 
Next consider the subregion B(i) in which 

a/r < 1 — y± < /i and a/r < y 2 — 1 < 1 — y\ (A3) 
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Expanding the integrand in (Al) given these conditions we have 



3/2 — 2/i 



1 



\(2/2 — 1)2/2(1 — 2/1)2/1/ (l-£/i) r V 
The corresponding contribution to the double integral is 



ri^i + ... 
2/2-1 



1— <r/r 



dyi 



l+(l-2/l) 



) r Jl+a/r 



dy 2 + T 



l—a/r 



dyi 



l-ix 



* r (l-2/l) r 

When evaluated for large-r these two integrals give 



l-2/i) 



r-i 



dy 2 



l+^/r 2/2-1 



(A4) 



r-2 



and 



(r-2) 2 [W 



r-2 



+ 



r-2 



r-2 



logr 



(A5a) 



(A56) 



respectively. The corresponding contribution to from (A5a) is 0(1/ (T — 2) 2 ) while 
the contribution from (A5b) is 0(1/ (T — 2) 3 ). Hence (A5a) does not contribute to C^\r). 

The analysis for region B(ii) proceeds analogously to that of region B(i) just presented. 
We again find contributions (A5a) and (A5b). In regions C(i) and C(ii) we find a con- 
tribution to A( 4 ) which is 0(1/ {T - 2) 2 ). This region therefore does not contribute to 
C^\r), and since all regions are now exhausted, the terms in the asymptotic expansion of 
p++^(r) which contribute to A*- 4 - 1 give a contribution to the latter which is 0(1/ (F — 2) 3 ). 
Adding all such terms, including the term corresponding to (A2) which gives a divergent 
contribution to A^, we obtain the expansion 



P++ (r) ~ C 
2r / 



1 



+ 



,2r-2 



(r-2) 



(T-2Y 



2 r 2F-2 



(1 



r \ 
a J 



r-2 , 



) 



r-2 



r-2 



logr 



(A6a) 



The first term originates from subregion A while the remaining two terms have equal 
contribution from subregions B(i) and B(ii). 

A similar approach suffices to analyse the terms in the asymptotic expansion of (2.23c) 
and (2.23b) which contribute to A^. We find 



P T - { :\r) 



C 4 r 



(r - 2) 2 r 2r ~ 2 



and 



P?- 4) (r) ~ C 



fr 



2)2 r 2r-2 



r-2 , 



r-2 



- 1 



(A66) 



+ 



sr 

,2F-2 



(v-2y 



r-2 



-11 



r-2 



r-2 



logr 



(A6c) 



The particle configurations corresponding to the subregions of integration contributing to 
the terms in these expansions are given in Eig. 3. 
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Forming the combination (2.24) of (A6a),(A6b) and (A6c) required to form C(r) we see 
that a lot of cancellation takes place, and the result (2.27) of the text results. Furthermore, 
it is important to observe that due to the cancellations only p+^(r) contributes to 
and this contribution is restricted to the region of integration depicted in b(ii) of Fig. 3. 
Analogous to (A5b) this contribution can be written as a double integral which is given 
by (2.28) of the text. The integrand of the double integral is obtained from the expansion 
of the integrand (excluding the factor of l/r r ) in the second double integral of (2.23b): 

/ x(y — r) \ r /l\ r r 



(y — x)y(x — r) J \y — x J (x — r) (y — x) 



r-i 



valid for x — r and y — x small and positive with y — x < x — r. The domain of integration 
is that implied in configuration b(ii) of Fig. 3. 



APPENDIX B 



Our objective here is to study the leading large-r behaviour of the formula (3.4) for C^(r) 
as a function of T. We begin by changing variables x\ wr«X2 ^ rv, so that (3.4) reads 



«V 4 



dudv 



u — v 



(v — l)(u — l)uv 



2S 



u(v — 1) 



(u — v)(u — l)v 



\u — v\ 



(Bl) 



where * denotes the integration region 

{(u,v) G R 2 : \u — v\ > cr/r, \u\, \v\ > a/r, \u — l\,\v — 1\ > a/r} 

Due to the integrand being symmetric in u and v, and unchanged by the mappings 
u \— > 1 — u, v i — ^ 1 — we can restrict the integration region to 

{(u,v) G R 2 : u > v, v > 1 — u, \u — v\ > a/r, \u — l\,\v — 1 1 > cr/r} (B2) 

provided we multiply the integral by 4. Let us divide the integration region (B2) into 
labelled and unlabelled regions as indicated in Figure 4. 

All regions outside those labelled in Fig. 4 give a contribution to C^(r) which is 
0(l/r 2T ~ 2 ). This follows from (Bl) since the range of the integrand is independent of r in 
these regions. The range of the integrand is not independent of r in the labelled regions 
of Fig. 4, and consequently further calculations are needed. 

Let us illustrate our method of analysing the leading order contribution from the 
labelled regions by considering region D. The total contribution to C^(r) from this 
region is given by 



g 2 C 4 

r,2r-2 



roo ( rn l—<r/r\ 

/ du + dvF(u,v) 

Jl+H \Ja/r J-n I 



where F(u, v) is the integrand in (Bl). For large-r the leading contribution to (B3) comes 
from the neighbourhood of \v\ = a/r. We can thus ignore the last term in (Bl), and 
the term implicit in the symmetrized form of the second term in (Bl). Expanding the 
remaining two terms for small-v gives 

,2 \ 



2r 



(u-iy 



v 1 f v\ 2 -ir 

■Z + V -2\Z) + T + °("). 
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in place of F(u,v) in (B3). 

Evaluating this integral to leading order for large-r gives behaviours 

°(^) 2 <r<3, o(^)r = 3, o( (r _3 )rr+1 ) T>3 (B4) 

Proceeding similarly, we find the leading large-r contributions from the labelled regions 
of Fig. 4 to be as in the table below. 



Region Contribution to C^(r) 

2<r<3 r=3 r>3 

A,B,C,D O(^) O(^) 0( jr ^ rTT ) 
O(^) O(^) 0(£) 

P o(^) o(^) o(^) 

Hence the leading order contribution to C^ 4 ^(r) comes from region a for 2 < T < 4 
and from region (3 for Y > 4. These regions correspond to the charge configurations of 
Fig. 5. 
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Figure 1. The two distinct nested dipoles at 0(( 6 ). The square braces indicate a 
screening operator which connects the charge on the left end with the dipole on the right 
end. 

Figure 2. Some regions of integration in the integral (Al) for p++^(r). 

Figure 3. The configuration (a) gives the term in (A6b), while the configuration b(i) 
gives the first term in (A6c). The configurations b(ii) and b(iii) both contribute equally 
to the remaining terms in (A6c). These configurations each have a multiplicity factor of 
4, which corresponds to the four ways of placing the mobile pair x and y on either side of 
the two root charges. 

Figure 4. Regions of integration in (B2) which are used to analyse the integral (Bl) for 
C^\r). A distance a/r either side of each dotted line bordering or contained within a 
region is to be excluded. 

Figure 5. Charge configurations corresponding to regions a and (3 in Figure 4, which 
give the leading order contribution to C^\r) for 2 < T < 4 and T > 4 respectively. The 
arrows indicate alternative placement of mobile charges. 
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